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1,0  INTRODUCTION 

This  report  summarizes  the  results  of  the  Doppler  Navigation  Study  conducted 
by  the  Submarine  Signal  Division  of  Raytheon  Company  for  the  Office  of  Naval  Re¬ 
search  under  Contract  NOnr  3376(00).  The  purpose  of  the  study  was  to  establish  the 
feasibility  of  applying  Doppler  principle  to  sonar  navigation  in  shallow  depths.  The 
study  encompassed  both  analytic  and  experimental  work,  including  the  fabrication  and 
testing  of  breadboard  circuitry  necessary  for  the  experimental  phase .  The  results  of 
the  study  prove  the  feasibility  of  shallow  water  sonar  Doppler  navigation. 

A  system  configuration  of  minimal  navigator  error  is  described  herein,  together 
with  the  results  of  experimental  error  studies.  This  study  has  generated  the  necessary 
engineering  information  for  design  of  a  shallow  water  sonar  Doppler  navigation  system. 
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2.0  DISCUSSION  OF  PRINCIPLES  AND  EXPECTED  PERFORMANCE 

This  section  describes  the  Doppler  effect  briefly,  the  geometry  of  the  problem, 
and  other  parameters  of  the  ejqierimental  apparatus  used  to  validate  the  analytic  re¬ 
sults  obtained  during  this  study. 


2.1  Doppler  Effect 


The  Doppler  effect  between  a  moving  source/observer  and  a  fixed  reflector  is 


accurately  described  by  the  following  equation; 

2 

.  .  V  V  v 

2  V  f  cos  a  /  L  w 


1  - 


?■] 


(2.1.1) 


V  Doppler  shift,  cps 

V  Ship*  s  horizontal  velocity,  fps 
f  =  Transmitted  frequency,  cps 

6  "  Angle  between  the  acoustic  axis  and  the  reflecting  surface 

(depression  angle) 

C  >  Velocity  of  sound  in  the  medium,  fps 
=  Velocity  of  the  medium,  fps. 

Because  of  viscous  friction  on  most  bottoms,  tiie  velocity  of  the  xnedium  at  the  reflecting 
surface  is  usually  small  and  thus  may  be  neglected  in  the  limit,  giving 

V  =  2Vf  cos  a/(U- V),  cps  (2.1.2) 

For  small  values  of  V,  relative  to  C,  C  -  V  %  C  giving  the  simplified  Doppler  equatim 
2Vf  -  “ 

V  =  -g-  cos  6,  cps  (2.1.3) 


Equation  (2.1.3)  gives  values  for  v  within  0.5%  for  V  less  than  14.3  knots,  as  com¬ 
pared  with  Equation  (2.1.2). 
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2.2  Balanced  Acoustic  Array 

In  the  general  sense,  the  Doppler  equation,  (2.1.3),  is  sensitive  to  any  motion 
relative  to  the  reflecting  surface,  not  only  horizontal,  but  also  vertical  motion  due  to 
bottom  slope  or  ship  heave.  In  order  to  remove  the  effects  of  vertical  components,  a 
balanced  or  'Janus'  system  is  used.  This  system  is  discussed  below. 

If  two  acoustic  receive-transmit  axes  (hereafter  called  acoustic  axes)  are  arranged 
in  a  vertical  plane  as  shown  in  Figure  1,  it  is  easy  to  show  that,  as  long  as  the  condition 
of  verticality  is  maintained,  the  difference  of  the  Doppler  shifts  of  the  acoustic  axes  is 
insensitive  to  vertical  components  of  motion. 


U 


The  shift  caused  by  the  combined  effects  of  horizontal  velocity,  V,  and  vertical 


velocity,  u,  along  axis  is 


2f 


vj  =  (V  cos  -  u  sin  dj)  — 


(2.2.1) 


and  the  shift  along  N  is 


v  =  -  (V  cos  d  +  u  sin  d  ) 

«j  o  3 


2f 


(2.2.2) 


Substracting  (2.2.2)  from  (2.2.1)  with  the  assumption  d  =  d  gives 
4V^fcosdg  ‘ 


1 


‘’3‘ 


(2.2.3) 
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Equation  (2.2.3)  indicates  that  the  Doppler  shift  obtained  by  a  balanced  system  is  insensi¬ 
tive  to  vertical  motion  and  is  twice  the  value  of  that  obtained  by  horizontal  motion  of  a 
single  axis.  It  may  be  said  without  formal  demonstration  that  these  advantages  accrue 
only  if  the  depressicm  angles  are  equal  with  respect  to  the  true  horizontal,  indicating 
that  best  accuracy  requires  that  the  acoustic  axes  be  space  stabilized. 

2.3  Error-Producing  Factors 

Errors  accrued  in  the  velocity  measuring  process  can  generally  be  broken  down 
into  the  categories  of  1)  Geometry,  2)  Medium  and  3)  Direction.  The  general  geometric 
parameters  can  be  well-defined,  as  can  the  paranoeters  involved  in  determining  directicm, 
but  there  are  certain  characteristics  of  the  medium  which  are  not  too  well-defined  at 
these  frequencies.^ 

2.3.1  Vertical  Standoff  Errors 

The  errors  due  to  a  vertical  standoff  and/or  pitch  and  roll  can  be  quite  easily  de¬ 
fined  if  certain  assumptions  are  made.  The  following  assumptions  will  be  made  in  this 
section;  1)  No  vertical  motion  effects  will  be  considered,  2)  The  depressim  an^es  in 
each  axis  are  the  same,  3)  The  separation  angles  between  transducers  are  exactly  90*, 

4)  All  time -variable  effects  are  symmetrical  and  occur  centrally  about  the  ideal  static 
situation  and  S)  All  static  offsets  occur  in  the  absence  of  c^amic  variations. 

Ccmsider  first  the  static  situation  where  the  array  vertical  has  a  standoff  error  witii 
respect  to  the  true  vertical,  the  standoff  angle  being  in  the  vertical  plane  containii^  the 
depression  angles.  If  the  trim  angle  is  6  as  shown  in  Figure  2,  the  Doppler  shifts  in  the 
1  and  3  axes  become 

V.  K  cos  (d -i-O)  -  K  cos  d(cos6 -tan  d  sinO)  (2.3.1) 


References  will  appear  on  page  35. 
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FIGURE  2  TRIM  ERROR  GEOMETRY 


where 


K  = 


-K  cos  (6  -6)  = 

2Vf 

C 


-K  COS  6  (cos  64^  tan  6  sin  6) 


The  error  is  given  as 

E,  =  100  (cos  e  -  tan  «  sine  -  i)% 


for  the  1  axis  and 

E^  ^  100  (cos  e  -  (1  +  sin  6)  ]  % 

for  6  -  45°  . 

For  the  balanced  system,  the  error  is 

E,,  =  -  100  (1  -  co8e)%  %  -  . 


which  is  independent  of  the  depression  angle. 


(2.3.2) 


(2.3.3) 


(2.3.3a) 


(2.3.4) 


Error  values,  plotted  for  both  the  single  and  Janus  cases,  are  shown  in  Figure  3. 
Note  that  the  error  is  always  negative  in  the  Janus  configuration,  due  to  the  asymmetry 
of  the  error  of  the  single  axis.  At  ±3°  inclination  the  Janus  error  is  -0.137%,  where  the 
single  error  is  5.09'7  for  a  negative  inclination  of  3°  and  -5.37%  for  a  positive  inclination 
of  3 '.  At  the  7  °  inclination  points  the  single  system  asymmetry  is  quite  evident,  the 
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VERTICALITY  ERROR  DEGREES 

FIGURE  3  ERROR  INTRODUCED  BY  NON-VERTICALITY  IN  JANUS  AND 
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difference  in  magnitude  of  the  errors  caused  by  the  positive  and  negative  inclination 
being  1.5%.  The  error  in  the  Janus  configuration  is  about  -0.75%  for  this  7*  inclina¬ 
tion.  The  superiority  of  the  balanced  arrangement  in  reducing  indicated  errors  caused 
by  array  inclination  is  quite  evident. 


2.3.2  Pitch  and  Roll  Errors 

The  effect  of  a  symmetric  motion  of  the  array  may  be  considered  by  allowing  the 

angle  6  to  become  a  function  of  time,  say  0(t)  =  6  sin  u  t  where  0  is  the  peak  devia- 

m  m 

tion,  or  amplitude,  of  the  motion  in  the  vertical  plan  and  u  is  the  angular  frequency  of 
the  motion.  For  the  single-axis  case  the  indicated  shift  may  be  integrated  in  time  over 
a  full  cycle  of  the  error  angle  6,  and  the  mean  value  of  the  indicated  shift  with  the  error 
angle  may  then  be  compared  with  the  indicated  shift  when  the  error  angle,  6,is  zero. 


If  the  single-axis  Doppler  shift  given  by  Equation  (2.3.1)  is  modified  ford  =  45* 
and  6  «  6(t), 

. .  [  cos  (6  Bin<>>t)-sin(6  sinwt)  ]  ,  (2.3.5) 

1  «  xn  m 

2  TT 

This  may  be  integrated  over  a  full  cycle  of  6  ,  which  corresponds  to  a  period  of  —  . 

m  w 

The  average  shift  over  this  interval  is  then 

T 


[cos  (6  sinwt) -sin  (6  sinwt)]dt 
m  m 


which  is 


Comparing  this  with  the  shift  obtained  when  0  =  o  -  6, 

m 

the  following  is  obtained  for  the  percentage  error; 


2/2  - 
K  ’'l 


Jo  (6) 


(2.3.6) 


Error  «  -100  [  I  -  J  (0  )  1  %. 

o  m 


(2.3.7) 
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The  same  sort  of  calculation  may  be  carried  out  to  llliutrate  the  effects  of  pitch  <ni 

a  balanced  system.  If  equation  (2.3.2)  Is  subtracted  from  equation  (2.3.1). 

v,  -  V,  =  2K  cos  6  cos  6 
X  3 

and  the  steady  shift  Is 


Vj  ~  Vg  =  2 K COS 6  ,6  =  0. 

The  average  over  a  full  cycle  of  6(t)  is  given  l^ 

IT 

W 

giving  the  same  error  for  sinusoidal  pitch  as  obtained  for  the  single  system, 

Error  =  -100  [  1  -  J  (0  )]  %  . 

o  m 


(2.3.8) 


(2.3.9) 


(2.3.7) 


The  error  values  given  by  equatiim  (2.3.7)  are  shown  in  Figure  4.  On  the  basis  of 
the  tabulated  values  in  Figure  4  it  can  be  seen  that  some  sort  of  compensation  or  stabili¬ 
zation  is  necessary  if  hi^ly  accurate  navigation  information  is  to  be  maintained  for  all 
weather  conditions.  It  mif^t  be  noted  that  the  error  produced  by  this  condition  always 
reduces  the  indicated  velocity. 


2.3.3  Sound  Velocity  Errors 

Equation  (2.1.3)  indicates  that  the  indicated  velocity  (Doppler  shift)  is  dependent  on 
the  velocity  of  sound  in  the  medium,  C.  If  this  velocity  is  not  known  and  properly  com¬ 
pensated,  there  will  be  a  seating  error  in  the  indicated  velocity.  For  practical  (Mirposes, 
this  error  appears  as  a  systematic  error,  probably  nearly  constant  over  a  period  of  many 
days  for  any  sheltered  location.  Since  the  major  cause  of  change  in  the  velocity  of  sound 
is  due  to  temperature  variations,  the  error  due  to  variation  in  temperature  has  been  com¬ 
puted  and  is  illustrated  along  with  the  corresponding  depression  angles  needed  to  keep  the 
shift  factor  constant  for  a  system  with  a  nominal  45*  depression  angle. 
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ANGULAR  ERROR  9,  DEGREES  ARC 


FIGURE  4  ERROR  IN  INDICATED  DISTANCE  TRAVELLED  DUE  TO  ARRAY  VERTICAL  OFFSET 
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The  velocity  of  sound  may  be  t^jproximated  Kuwahara' s  formula ; 


where 


C  =  4422  +  11.25  T  -  0.045  T  +  4.3(8-34) 

C  =  Velocity  of  sound,  fps 
T  =  Temperature,  “F 
S  =  Salinity,  parts  per  thousand 
and  the  simplified  Doppler  equation  may  be  written  as: 


(2.3.10) 


2  Vgf  cos  6 


V  - 


4422  +  11.25T  -0.045T 


(2.3.11) 


if  8  is  approximately  34  parts/1000  and  the  value  of  depth  is  small.  If  v/V  is  fixed  at 

s 

some  convenient  value,  say  100  cps/knot,  then  the  desired  depression  angle,  6,  may  be 
computed,  knowing  f  and  T,  as 


6  =  cos 


-1 


2Vgf 


(4422  +  11.25T  -  0.045  T  ) 


(2.3.12) 


Values  of  5  computed  for  a  nominal  6  of  45*,  f  of  210  kc  and  v/V  of  100  cps/knot  = 
59.25  cps/fps  are  shown  in  the  upper  portion  of  Figure  5.  The  lower  curves  of  Fig¬ 
ure  5  illustrate  the  velocity  C  versus  water  temperature  as  computed  from  equation 
(2.3.8).  Computed  values  of  indicated  error  are  for  an  uncompensated  situation. 
Fifty  degrees  Fahrenheit  was  chosen  as  a  mid-point  and  the  error  was  computed  for 
every  10"F  from  30*  to  80*F  using 


100 


^50  -  ^x 
^50 


=  Error 


(2.3.13) 


2.3.4  Linearity 

As  mentioned  in  8ection  2.1  the  linear  Doppler  equation,  equation  (2.1.3),  has  an 
error  of  approximation  associated  with  it  when  compared  with  equation  (2.1.2).  The 
magnitude  of  this  error  may  be  easily  determined  using  the  relatimi 

1 


Velocity  Error  =  lOOjl  • 


% 


(2.3.14) 


1  - 
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where  V  may  assume  either  a  positive  or  negative  value  according  to  whether  V  is  in 
the  direction  of  propagation  (positive  V)  or  against  the  direction  of  propagation  (nega¬ 
tive  V).  Using  an  assumed  velocity  of  sound  of  4800  ft-'sec,  the  velocity  error  for  a 
single  axis  system  is  tabulated  below.  It  can  be  seen  that  this  error  is  quite  large  for 


the  hi^er  ship  velocities. 


TABLE  I  -  VELOCITY  ERROR 


Speed,  Single  System  Balanced  System 

Knots  V^elocity  Error  Velocity  Error 


1  0.035% 
10  0.352% 
50  1.76% 


1.24(10“®)% 

0.00124% 

0.031% 


For  the  balanced  system,  the  error  due  to  this  factor  is  small,  the  error  being 
given  by; 


Balanced  Velocity  Error 


(2.3.15) 


values  for  which  are  tabulated  in  Table  I  also.  These  figures  indicate  that  a  balanced 
system  is  definitely  advantageous  on  hi^er  speed  ships  if  no  compensation  for  the  velo 
city  error  is  to  be  computed  by  the  Doppler  navigation  device. 


2,3.5  Compass  Error 

The  primary  error  source  in  a  stabilized,  balanced  Doppler  Navigation  system  is 
compass  error,  or  heading  reference  error.  The  magnitude  of  error  introduced  by  the 
compass  can  be  expressed  in  either  the  cross-track  and  track  directions  or  in  earth 
coordinates.  Figure  6  indicates  both  errors. 

The  true  coordinates  as  represented  on  the  earth  and  on  the  navigator  displays  are 
the  ordinate  and  abscissa  of  Figure  6.  Line  OA  is  the  actual  track,  or  velocity  vector  of 
length  A,  The  coordinate  system  N^OE^  is  the  rotated  coordinate  system  produced  by  the 
compass  when  an  error  angle  E^  exists.  Knowing  that  OA  is  the  same  length  in  either  co¬ 
ordinate  system,  the  errors  may  be  quite  simply  derived  from  the  compass  error,  E  ,  and 

c 

the  course  made  good,  C. 
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The  northerly  error  is 

E„  =  ON'  -  ON  =  A  cos  (C  -  E 
N  c 


)  -  A  cos  C 


and  similarly  the  easterly  error  is 


(2.3.16) 


E_  =  OE'  -  OE  =  A  sin  (C  -  E  )  -  A  sin  C  (2.3.17) 

£  C 

The  normalized  coordinate  errors  are  not  particularly  useful  since  one  or  the  other  will 
be  undefined  if  C  is  a  multiple  of  t/2. 


The  cross-track  error,  E^.j„  and  the  track  error,  E.^,,  are  more  useful  in  the  norma¬ 
lized  form,  since  they  do  not  depend  on  the  course,  C.  Since  the  indicated  track  OA'  and 

the  actual  track  OA  form  an  isoceles  triangle  A'  OA  with  vertex  angle  E  ,  the  line  A'  A  is 

c 

the  chord  of  a  circle  of  radius  A  defined  by  a  central  angle  E  .  This  chord,  E,  is 

E  ® 

c 

E  =  2  A  sin 

in  length  and  the  triangle  A'  AB  may  be  solved  to  give  cross-track  error 

E  EE 

E  =  E  cos  =  2  A  sin  cos  (2.3.18a) 

X  1  2  it  it 


E_  =  AsinE  5;AE 
xT _ c  c 

The  error  along  the  track  direction,  E.^ ,  is  given  by 

e  2  e 

E_  =  -E  sin  ~  =  -2  A  sin  — 

T  2  2 


(2.3.18b) 


(2.3.19a) 


or 


-  A  (1  -  cos  e)  % 


(2.3.19b) 


which  is  much  smaller  than  E  _ .  Since  E  _  is  much  larger  than  E_  and  since  the  other 

xT  xT  o  -j. 

error  factors  tend  to  act  along  the  track  if  the  direction  of  travel  is  close  to  a  cardinal 
course ,  the  compass  error  may  be  estimated  if  a  traverse  between  two  known  points  on  a 
cardinal  line  is  made.  This  fact  has  been  utilized  in  calibrating  the  orientationof  the  experi¬ 
mental  equipment  in  lieu  of  trying  to  boresi^t  the  transducers  to  the  cardinal  points,  a 
task  which  could  become  quite  complicated. 
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2.3.6  Spectral  Contributions 


The  effect  of  finite  aperture  on  the  spectrum  of  back-scattered  Doppler  information 
has  been  the  subject  of  many  articles  in  the  literature  of  Radar  Doppler  navigation.  The 
bibliography  in  the  last  section  of  this  report  lists  some  pertinent  documents  on  the  gen¬ 
eral  topic  of  Doppler  navigation. 

3 

Berger  (page  109)  gives  the  width  of  the  Doppler  spectrum  as 


2  V 


(sin  a  )  A  6 
o  o 


where  A6  =  round  trip  half -power  width  of  the  beam  pattern, 
this  study,  allowing  for  the  Janus  scale  factor,  this  gives 
Av  =  200V  A  a  ,  V  in  knots 


For  the  system  used  in 


The  half -power  point  of  both  receiving  and  transmitting  units  of  the  experimental  system 
is  about  6* ,  giving  an  effective  round  trip  width  of  4. 25*  =  0.074  radians,  leading  to  an 
expected  value  of 


^  149  ^ 

V  ^  ^  '  knot 


or  74  cps  for  5  knots  velocity. 


A  series  of  spectral  measurements  were  taken  with  the  transducers  locked  in  the 
forward  direction.  The  forward  speed  was  held  at  5  knots  nominal  and  water  depth 
varied  from  20  to  160  feet.  The  analyzer  used  was  the  multiple -filter  Rayspan  with 
filter  width  of  8  cps  and  the  sampling  rate  was  200  per  second.  The  spectral  width 
varied  from  50  cps  to  75  cps;  the  spectral  width  was  estimated  from  20  scans  of  the 
filter  bank  photographed  on  the  face  of  an  oscilloscope.  Unfortunately,  tlw  camera  has 
a  li^t  leak,  rendering  most  of  the  photographs  unsuitable  for  reproduction;  Figure  7  is 
a  tracing  of  the  envelopes  of  spectra  taken  in  30  and  150  feet  of  water.  Continuous  obser¬ 
vation  of  the  spectra  while  traversing  the  range  of  depths  did  not  indicate  any  apparent 
correlation  between  spectrum  width  and  depth.  Speed  changes  and  rapid  depth  changes 
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did  appear  to  widen  the  spectrum  slightly  but  by  no  more  than  a  factor  of  1 . 5.  This  is 
a  difficult  point  to  interpret  since  the  whole  spectrum  slides  with  speed  changes,  and  depth 
changes  would  be  expected  to  act  differently  on  the  fore  and  aft  beams  over  the  short  sampl¬ 
ing  interval  (3  milliseconds). 

The  observed  spectra  correspond  very  closely  to  the  expected  value  of  74  cps  for 

a  velocity  of  5  knots  as  computed  previously.  The  expected  contributicm  of  spectrum 

3 

width  to  system  error  may  be  computed  fay  the  method  given  by  Berger  (pages  110,  111) 
and  others. 

Using  Berger*  s  results: 

H.  -  ""P  .V,  _L  =  i- 

V  y  2T  2  cos  6  "W' ' 

=  normalized  RMS  velocity  error, 

<rD 

~  =  normalized  RMS  distance  error, 

T  =  ^  =  time  to  traverse  distance  D, 

an  expected  value  of  o-  D/D  for  D  ■  1  mile  may  be  conq>uted  using  the  following  para¬ 
meters; 

6  =  45* 

/  =  c/f  =  0.02285  ft, 

D  =  6076.1  ft. 

Ad  =  0.074  radians, 

<r  D 

giving  =  0.03%.  The  normalized  value  vD/D  varies  inversely  with  the  square  root  of 
distance  so  that  it  can  be  seen  that  the  expected  error  c<mtribution  of  the  spectral  distribu¬ 
tion  of  the  Doi^ler  signal  soon  becomes  insignificant  when  compared  with  other  error  fac¬ 
tors,  such  as  the  gyrocompass,  which  can  be  expected  to  contribute  errors  of  the  order 
of  0. 1%  due  to  composite  errors  in  the  data  transmission  system  alone! 
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4 

Note;  Fried  ,  (Appendix  I)  derives  an  expected  Janus  spectrum  width, 
based  on  Berger' s  method,  which  is  one -half  the  width  of  that 
given  here.  A  performance  factor  is  included  which  in  later  sec¬ 
tions  of  Appendix  I  of  this  reference,  is  taken  as  2,  thus  result¬ 
ing  in  the  equivalence  of  the  two  systems  when  considered  in 
practice . 

There  have  been  predictions  of  observed  errors  of  a  larger  magnitude  than  the  above 
due  to  the  shift  in  the  apparent  power  center  of  the  ensonified  region  of  the  bottom,  due 
to  the  change  in  the  exponent  of  the  scattering  law  as  regions  of  differing  acoustic  im- 
pedance  are  traversed  .  Since  no  measurements  of  backscattering  have  been  taken  in 
this  study,  no  conclusive  data  may  be  offered,  but  operation  in  both  Narragansett  and 
Chesapeake  Bays  have  not  produced  errors  which  would  seem  to  be  ascribed  to  this 
source. 

2.3.7  Conclusions  Regarding  Error 

The  major  error  factors  acting  on  a  sonar  Doppler  navigator  are  the  heading  and 
sound  velocity  errors.  Of  these  two,  the  heading  error  is  the  most  serious  since  the 
control  of  heading  reference  error  is  the  most  difficult  to  predict  and/or  correct.  It 
is  anticipated  that  the  most  effective  method  of  controlling  the  sound  velocity  error  will 
be  acconq)lished  through  the  use  of  a  sound  velocity  controlled  carrier  frequency  (con¬ 
stant  wavelength  transmission).  To  accomplish  this,  the  tuned  sense  discriminators 
must  be  replaced  by  a  discriminator  that  is  not  dependent  on  a  fixed  transmission  fre¬ 
quency. 

It  has  been  noted  experimentally  that  the  relationship  of  equation  (2.3.10)  between 
sound  velocity  and  depression  angle  is  not  exact,  probably  due  to  the  fact  that  the  energy 
is  distributed  in  a  finite  beam  width  rather  than  in  the  infinitesimal  beam  assumed  in  the 
elementary  analysis  of  the  problem.  Due  to  this  fact,  and  considering  the  mechanical 
tolerances  that  would  have  to  be  held  in  compensating  the  sound  velocity  by  changing  the 
depression  angle ,  this  method  of  compensation  aiq)ears  to  be  impractical  for  anything  but 
a  laboratory  device .  The  aforementioned  constant  wavelength  transmission  method  would 
be  recommended  for  an  operatimial  navigator . 
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The  analysis  outlined  in  the  Second  Interim  Report  on  this  project  was  completed 
and  compared  with  experimental  results,  but  in  all  instances,  the  actual  performance 
was  better  than  predicted. 

This  is  undoubtedly  due  to  the  fact  that  continuous  sinusoidal  pitch  and  roll  motions 
are  an  extreme  worst  case  condition.  The  use  of  the  data  of  Figure  4  gives  adequate 
prediction  of  the  errors  to  be  expected  due  to  imperfect  stabilization  when  the  average 
roll  amplitude  is  taken  as  the  independent  variable .  This  is  a  much  more  economical 
method  of  prediction  than  solving  the  complete  three -axis  problem  for  each  set  of  con¬ 
ditions  involved. 

In  any  case,  the  mean  square  error  contributed  by  the  gyrocompass  will  predomi¬ 
nate  in  a  stabilized,  compass  and  sound  velocity -compensated  navigator.  It  may  be 
stated  then  that  for  all  practical  purposes,  shallow  water  Doppler  navigation  accuracy 
is  at  present  limited  by  the  accuracy  of  available  heading  references. 


CONFIDENTIAL 


19 


R332 


CONFIDENTIAL 


3.0  EXPERIMENTAL  RESULTS 

This  section  describes  the  results  obtained  with  three  different  configurations  of 
equipment.  The  errors  of  closure  are  presented  for  ninety-five  runs,  of  which  the 
majority  were  2.7  miles  in  length.  Scaling  errors  are  not  considered  because  of  the 
lack  of  exact  data  on  shore  positions  and  the  fact  that  sound-velocity  calibration  and 
compensation  were  not  accomplished  on  a  continuous  basis.  Demonstration  of  closure 
capability  indicates  the  attainable  performance  when  scaling  parameters  are  properly 
adjusted.  The  operating  depths  on  these  cases  ranged  from  2  to  30  fathoms. 

3.1  Case  I  -  Unbalanced  System 

Fifty-eight  runs  were  made  using  an  unbalanced,  unstabilized  navigator.  No  com¬ 
pass  compensation  was  applied.  This  corresponds  to  what  would  now  be  considered  a 
sub-marginal  operational  system.  The  following  histograms  show  the  distribution  of 
errors  (Figure  8).  The  statistics  computed  for  this  series  are; 

Mean  North-South  Error  0.13%  North 

Standard  Deviation,  North-South  0.95% 

Mean  East -West  Error  0.6%  East 

Standard  Deviation,  East -West  3.6% 

The  existence  of  the  distinct  correlation  between  North  and  East  errors  may  be  attributed 
to  the  gyrocompass  since  the  starting  transient  was  to  the  northwest,  which  would  be  ex¬ 
pected  to  produce  a  westerly  compass  error  and  an  easterly  indicated  error  on  the  leg  of 
the  run  which  was  predominately  North-South,  and  a  northerly  error  on  the  leg  which  was 
predominately  West-East.  Althouj^  it  was  not  detected  at  this  time,  later  tests  showed 
that  rough  water  caused  a  severe  distortion  of  Doppler  signal  when  running  at  speeds 
above  5  knots.  This  effect  was  correlated  with  the  appearance  of  large  numbers  of  air 
bubbles  in  the  sea  chest  of  the  training  gear,  thus  pointing  toward  the  conclusion  that 


CONFIDENTIAL 


20 


R332 


CONFIDENTIAL 


bubbles  from  the  bow  wave  were  being  introduced  into  the  acoustic  path  of  the  navigator. 
This  may  have  been  a  contributing  factor  toward  the  extreme  peak  errors  noted  in  Cases  I, 
n  and  in,  but  since  this  effect  was  not  noticed  until  after  the  data  were  taken,  the  data 
were  not  changed. 

3.2  Case  II  -  Balanced  System 

The  acoustic  configuration  was  modified  to  the  Janus  form  (axes  in  all  cardinal 
directions)  to  eliminate  possible  errors  due  to  depth  variations.  The  Mark  19  gyro¬ 
compass  was  installed  without  velocity  correction  and  the  acoustic  and  electronic  equip¬ 
ment  was  installed  without  physical  stabilization  of  the  transducer  array. 

At  the  first  glance ,  the  results  were  rather  disappointing,  as  can  be  seen  in  the 
histogram  (Figure  9).  The  errors  were  not  as  large  on  the  whole  as  before  and  the 
variability  was  reduced,  but  the  design  goal  of  0.5%  was  not  attained.  A  consistency 
in  the  amount  and  direction  of  error  had  at  least  been  attained. 

The  statistics  for  this  case  were; 


Mean  error,  North-South 

1.52% 

Standard  Deviation 

1.28% 

Mean  error,  East-West 

1.05% 

Standard  Deviation 

0.47% 

Since  the  number  of  samples  (17)  was  small,  the  population  deviation  has  been  connq;>uted 
from  the  chi-square  distribution  at  the  95%  level  as 

0.183%  <_  o-NS  ^  0.765% 

0.067%  <  <rNW  <  0.28% 
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3.3  Case  ni  -  Velocity-Corrected  and  StablllMd  Balanced  System 

Velocity  correction  for  the  gyrocompass  was  Installed  in  May  1962  and  three  runs 
were  made  with  the  digital  display  disassembled,  thus  the  actual  errors  of  these  runs 
were  not  known,  but  the  charted  runs  showed  a  great  reduction  in  the  tendency  of  genera¬ 
tion  of  errors  only  northward  and  eastward.  The  stabilizer  for  the  transducer  array  was 
then  installed,  bringing  the  equtyment  to  its  present  state  and  the  errors  shown  in  Fig¬ 
ure  10  have  been  noted. 

The  statistics  for  this  case  were: 

Mean  error  -  North-South  0.078% 

Standard  Deviation  North -South  0.072% 

Mean  error  East-West  0.064% 

Standard  Deviation  East- West  0.096% 

As  for  Case  11,  the  population  deviation  has  been  computed,  on  the  basis  of  the  small 
sample,  to  be; 

0.0098%  <  <rNS  <  0.036% 

0.054%  <  «r  EW  <  0.198%. 

The  average  error  is  noticeably  less  and  the  variability  has  been  reduced  by  a  significant 
amount.  Certain  difficulties  axe  still  being  experienced  with  the  gyrocompass  in  that 
rapid  maxieuvering  produces  accelerations  which  result  in  heading  errors.  Extremely 
ti{^t  turns  at  7  or  8  knots  seem  invariably  to  cause  closure  errors. 
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4.0  EXPERIMENTAL  EQUIPMENT 

The  following  paragraphs  describe  the  equipment;  descriptions  of  some  parts  of 
the  equipment  were  presented  in  the  Third  Interim  Report  of  this  study. 

4.1  Acoustic  Equipment 

4.1.1  Transmitting  Transducers 

The  active  elements  of  the  transmitting  transducers  are  circular  disks  of  poled 
barium  titanate  ceramic  approximately  four  inches  in  diameter  and  approximately  one- 
half  inch  thick.  The  disks  operate  in  the  thickness  mode  of  a  frequency  of  210  kilo¬ 
cycles  per  second.  Being  over  14  wave  lengths  in  diameter,  the  directivity  index  is 
33  dB.  The  beam  width  at  the  3  dB  down  points  is  less  than  4. 5  degrees. 

Each  disk  is  power-factor  corrected  by  a  toroidal  core  transformer.  The  primary 
inqwdance  of  the  transformer  is  ai^roximately  32  ohms;  the  four  transmitting  trans¬ 
ducers  are  paralleled  to  obtain  an  8 -ohm  load  impedance  for  the  transmitter. 

The  transmitting  efficiency  of  these  units  is  in  excess  of  50  percent.  Each  unit 
receives  two  watts  of  electrical  power,  resulting  in  an  acoustic  ou^ut  power  in  excess 
of  one  watt  for  each  unit. 

The  transmitting  ceramic  disk  is  housed  in  a  cylindrical  aluminum  cup,  which  is 
anodized  for  corrosion  protection.  Acoustic  isolation  material  of  cork  in  a  neoprene 
binder  (corprene)  is  used  on  the  side  and  back  faces  of  the  disk.  The  disk  is  sig)ported 
in  the  housing  on  blocks  of  corprene;  the  housing  is  filled  with  an  elasticized  epozy  com¬ 
pound  to  provide  a  waterprocrf  assembly.  The  electrical  cable  is  led  throui^  a  threaded 
stud  which  serves  as  the  mounting  member,  and  is  located  on  the  bottom  d  the  aluminum 
cup. 

The  four  matching  transformers  for  the  transmitters  are  located  in  a  separate  alumi¬ 
num  housing  and  watertightness  is  assured  by  a  gasketed  cover  plate  and  silicone  ridiber 
encapsulaticm  of  the  transformers. 
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4.1.2  Receiving  Transducerg 

The  active  elements  of  the  receiving  transducers  are  identical  to  those  of  the  trans¬ 
mitting  transducers.  Unlike  the  transmitting  units,  however,  the  receiving  units  en4>loy 
a  novel  mounting  of  the  active  element. 

The  difference  in  transducer  construction  between  transmitting  and  receiving  imits 
arises  from  the  need  to  locate  the  receiving  preamplifiers  as  close  to  the  receiving 
transducer  as  possible.  For  this  reason,  the  preamplifiers  are  placed  in  the  transducer 
housing.  The  housing  is  basically  the  same  as  the  transmitter  housing,  being  an  anodized 
aluminum  cup.  Electrical  cable  is  led  throu{^  the  same  type  of  threaded  stud  on  the 
bottom  of  the  cup.  The  ceramic  disk  is  not  embedded  in  potting  compound  in  the  cup  as 
in  the  transmitting  units;  instead,  the  receiving  ceramic  disk  is  bonded  to  a  thin  anodized 
aluminum  dome .  The  periphery  of  the  aluminum  dome  fits  into  a  groove  on  the  upper  edge 
of  the  cup,  the  groove  then  being  filled  and  sealed  with  lead.  The  resulting  transducer 
assembly  leaves  sufficient  space  in  the  interior  of  the  cup  to  accommodate  the  preampli¬ 
fier.  The  preamplifier,  which  provider  45  dB  of  voltage  gain,  is  encapsulated  for  physi¬ 
cal  protection. 

4.2  Mechanical  Equipment 

4.2.1  Transducer  Array  Fixture 

The  transducer  array  fixture  has  four  purposes; 

a)  to  locate  the  four  transmitting-receiving  transducer  pairs  in  two  mutually 
perpendicular  vertical  planes. 

b)  to  locate  each  receiving  transducer  at  the  same  depression  angle. 

c)  to  provide  for  adjustment  of  the  depression  angle  for  temperature  com¬ 
pensation. 

d)  to  provide  maximum  acoustic  isolation  between  any  one  receiving  transducer 
and  all  other  transducers,  but  especially  the  associated  transmitting  trans¬ 
ducer. 
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Since  the  array  fixture  is  not  shielded  by  a  dome,  the  shape  of  the  fixture  must 
minimize  water  drag  forces.  The  shi^M  chosen  to  fulfill  all  of  these  requirements  is 
roughly  that  of  an  hourglass.  The  lower  inverted  cone  holds  the  four  transmitting  trans¬ 
ducers,  while  the  upper  inverted  cone  holds  the  four  receiving  transducers.  Joining  the 
two  Inverted  cones  is  an  upright  cone.  The  cones  position  the  transducers  while  the  place¬ 
ment  of  the  receiver  cone  over  the  transmitter  c(»e  provides  isolation  of  the  transmitter 
beam  from  the  receiver  beam. 

The  fixture  is  hollow;  each  transducer  is  located  inside  looking  out  through  a  circular 
hole  in  the  cone  face.  The  transducers  are  mounted  on  plates  pivoted  on  a  horizontal  axis, 
the  pivot  points  being  adjustable  so  that  the  transducers  can  be  aligned  to  look  in  mutually 
perpendicular  vertical  planes.  The  receivers  are  locked  together,  with  an  adjusting 
mechanism  to  provide  adjustment  of  the  depression  angle .  Adjustment  for  temperature 
compensation  is  affected  from  tq^iside  by  means  of  push-pull  flexible  cables  attached  to 
the  receiver  linkage.  As  the  cables  are  positioned  by  a  topside  screw  adjustment,  the 
linkagpe  produces  a  corresponding  rotation  of  the  receiving  transducers  about  their  hori¬ 
zontal  pivot  axis. 

Acoustic  isolation  between  transducers  is  effected,  first,  by  the  separaticm  of  trans¬ 
mitting  and  receiving  cones  as  mentioned  above  and,  second,  l^  g^enerous  use  of  isolation 
material  (Corprene)  to  break  any  acoustic  propagation  path  other  than  the  exterior  water 
path.  Each  transducer  is  covered  hy  corprene  1/4-inch  thick  on  the  transducer  case  to 
virtually  eliminate  transmission  or  reception  of  acoustic  energy  through  the  case.  Cor¬ 
prene  was  also  ^ipTied  to  the  exterior  face  of  the  upri^t  cmie  to  minimize  spuriotu  trans¬ 
missions  through  the  cone  into  the  receiver  beams. 

The  fixture  is  made  of  aluminum  plate  anodized  for  corrosion  protection.  At  various 
stages  in  the  progp'am  the  fixture  was  removed  from  the  water  for  modifications.  At  each 
removal  several  man-days  were  spent  removing  marine  growth  from  the  head  and 


CONFIDENTIAL 


27 


R332 


CONFIDENTIAL 


especially  the  transducer  faces.  Corrosion  of  the  head  became  noticeable  after  several 
months.  To  eliminate  these  problems,  an  antifouling,  anti-corrosion  protective  coating 
of  a  Raytheon  proprietary  material  was  applied.  After  a  six -month  immersion  following 
the  application  of  this  coating,  no  fouling  or  further  corrosion  was  noted. 

4.2.2  Azimuth  Drive  System 

The  transducer  array  fixture  is  rotated  by  the  azimuth  drive  system.  The  fixture 
is  fixed  to  a  vertical  shaft  protruding  throu^  a  sea  chest  in  the  hull.  (Interposed  be¬ 
tween  the  shaft  and  the  fixture  is  the  stabilization  system  described  in  the  next  section. ) 
The  shaft  rotation  is  controlled  l^  a  servomechanism  whose  reference  is  a  Mark  19  gyro¬ 
compass  .  The  shaft  is  rotated  so  as  to  maintain  the  north-south  transducer  axis  in  the 
north-south  direction,  as  indicated  by  the  g}rrocompass . 

The  control  or  error  s^cn&l  for  the  servomechanism  is  produced  by  comparing  the 
shaft  position  with  the  indicated  north  of  the  gyrocompass.  The  comparison  is  done  in 
a  two-speed  (1:1,  36:1)  synchro  transformer  mechanism  geared  to  the  shaft.  The 
error  signal  is  amplified  to  drive  a  motor  geared  to  the  training  shaft  to  null  the  error 
signal. 

4.2.3  Stabilization  System 

Interposed  between  the  training  shaft  and  the  transducer  array  fixture  is  a  stabili¬ 
zation  system.  The  purpose  of  the  stabilization  system  is  to  maintain  the  fixture' s  axis 
vertical. 

Basically  the  system  consists  of  two  horizontal  plates  Joined  by  a  flexible  vertical 
link,  which  is  a  ball  joint  inside  a  bellows.  This  construction  was  used  for  the  vertical 
link  to  achieve  a  constant  velocity  Joint  betwe«i  the  training  shaft  axis  and  the  transducer 
array  fixture  axis.  The  attitude  the  lower  plate  fixed  to  the  top  of  the  array  fixture  is 
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sensed  by  a  pair  of  electrolytic  levels  mounted  at  right  angles,  being  aligned  with  the 
North-South  and  East-West  axes  of  the  array  fixture.  A  deviation  from  the  vertical  pro¬ 
duces  an  electrical  signal  from  the  levels.  This  signal  is  amplified  and  applied  to  electro- 
hydraulic  servo  valves  which  then  apply  hydraulic  pressure  to  pairs  of  hydraulic  cylinders. 
Actuation  of  the  cylinders  produces  the  necessary  corrective  motion  to  bring  the  lower 
plate,  and  thus  the  array  fixture,  to  vertical.  The  servo  valves  and  the  hydraulic  power 
source  are  fixed  to  the  top  of  the  training  shaft. 

Although  the  purpose  of  the  stabilization  system  is  to  maintain  the  array  vertical 
independent  of  vessel  rolling  and  pitching,  the  stabilization  system  only  stabilizes  for 
small  vessel  motions.  For  larger  motions,  the  systemwillbut  attenuate  the  motions  of 
the  array.  The  fault  of  the  system  lies  with  the  electrolytic  levels  used  to  sense  devia¬ 
tions  from  vertical.  The  levels  are  subject  to  shifting  of  the  liquid  under  acceleration. 

The  resulting  system  tends  to  oscillate  at  a  frequency  in  the  region  of  the  frequency  vessel 
motion.  The  average  vertical  of  the  array  fixture  is  maintained,  however,  under  the 
li{^t  sea  conditions  normally  encountered  in  the  test  waters,  the  stabilization  system 
performing  satisfactorily. 

4.3  Display 

4.3.1  Plotter 

The  vessel' s  track  is  displayed  on  standard  navigational  charts  by  an  X-Y  plotter. 

The  plotter  employs  two  rods  at  ri^^t  angles  step-driven  by  rotary  solenoids  driving 
ball  nuto  along  lead  screws.  A  marking  pen  is  positioned  by  the  intersection  of  the  two 
rods.  Thus  as  the  rods  are  driven  in  the  X  and  Y  directions  in  correspondence  to  the 
input  pulses  in  those  directions,  the  pen  draws  the  track  of  the  vessel. 

Each  input  pulse  to  the  rotary  solenoid  produces  a  rod  step  motion  of  0.0073  inches, 
corresponding  to  a  travel  of  0.001  nautical  miles  on  a  chart  scale  of  1 ;  10, 000.  The 
direction  and  scale  factors  are  set  by  the  discriminator  and  scaling  circuits  as  indicated 
previously. 
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4.3.2  Countera 

The  distance  travelled  is  displayed  by  four  counters,  one  counter  for  each  cardinal 
direction.  These  counters  are  fed  the  same  pulses  as  the  plotter  but  are  mechanically 
independent.  The  counters  serve  to  record  the  total  distance  travelled  in  each  cardinal 
direction. 

4.3.3  Speed  Meters 

The  North-South  and  East-West  velocity  of  the  vessel  is  displayed  two  speed 
meters.  The  direction  of  the  vessel's  velocity,  whether  North  or  South,  is  indicated 
by  lights  activated  by  the  discriminator. 

4.3.4  Vessel  Heading 

The  vessel' s  heading  is  indicated  by  a  gyro  repeater  unit. 

4.4  Electronics 

The  electronic  portion  of  the  navigator  has  been  described  in  detail  in  the  Third 
Interim  Report  of  this  study  and  with  one  exception  remains  the  same. 

In  order  to  compensate  for  velocity  effects  on  the  gyrocompass,  a  compass  speed 
converter  unit  has  been  added  to  the  navigator. 

This  unit  converts  the  Doppler  frequencies,  North-South  and  East-West,  into  400- 
cycle  signals  scaled  to  provide  the  proper  control  to  the  compass.  The  Doppler  pulse 
information  drives  a  magnetic  pulse  conditioner,  the  output  of  which  is  rectified,  inte¬ 
grated  and  chopped  to  400-cycle  pulses.  These  are  amplified  and  filtered  to  provide  the 
proper  amplitude  and  waveform  for  the  conopass.  The  {Aase  of  the  400-cycle  signal  is 
controlled  by  switching  the  chopper  drive  {diase  with  the  direction  sense  signal  from  the 
discriminators.  The  scale  factor  has  been  set  at  0.75  volts  per  knot  which  is  fed  into  the 
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Mark  19  compass,  replacing  the  output  of  the  compass  speed  resolver.  Since  the  speed 
information  is  already  resolved  into  North  and  East  components  and  is  bottom -referenced, 
the  drift  setting  and  speed  resolver  of  the  compass  are  not  needed.  Since  the  speed 
deviation  angle  in  the  vicinity  of  Newport,  R.I. ,  amounts  to  0.08  degrees  per  North 
knot,  the  error  induced  by  uncompensated  speed  is  significant.  The  reduction  in  error 
between  the  second  and  third  series  of  runs  mentioned  before  is  thought  to  have  been 
effected,  in  good  part,  by  the  addition  of  this  speed  compensation. 

4.5  Significant  Equipment  Parameters  Summarized 

The  following  is  a  tabulation  of  the  significant  parameters  of  the  experimental 
equipment. 

4.5.1  Electrical 


Operating  Frequency 

210  Kilocycles/ second 

Transmission  Mode 

CW 

Electrical  Power  Output 

8  watts  at  210  kc 

Hydrophone  Preamplifiers 

Gain 

45  dB 

Bandwidth 

4  kc  to  clipper  inputs 

Clipper  Amplifiers 

Gain  to  clip 

112  dB 

Bandwidth 

80  kc 

Clip  Level 

-100  dB 

Output  Level 

2.5  V  RMS  (Saturated) 

Mixer  Bandwidth 

80  kc 

Low  Pass  Filter  Bandwidth 

1500  cpc 

Doppler  Voltage  Level 

2  VRMS 
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Scaling  Circuits 

Scale  Factors 


360,  1:5,000  | 

720,1:10,000  1 

1440,  1:20,000  [ 

2880,  1:40,000  j 

Discriminator 

Hysteresis 

2  cps 

Band  Limits 

±2  kc 

Con^iass  Speed  Converter 

Pulse  Input 

4  V 

Pulse  Input  Rate 

0  -  2  kc 

Output 

400  cps  sine  wave 

Output  Scale  Factor 

7.5  mV/cps 

Linearity 

<1% 

Null 

<100  mV 

Power  Iiq>uts 

115-V  60-cps  Single  phase 

!k750  W 

115-V  400-cps  Single  phase 

low  (choiqper  drive) 

440-V  60-cps  3-Phase 

1.5  H.  P.  motor  load 

Level  Control  An4>lifiers 

hiput 

60-cps  signal  from  levels 

Input  Level  Gradient 

%  1  Volt/degree 

Magamp  Output  Saturation  Level 

%  8  Volts  dc  into  1000  ohms 
servo  valve 

Gain 

10  (variable) 

Asimuth  Control  Amplifier 

Input  Signal  Source 

23CT6  sjrnchro  (1  ft  36  speed) 

Input  Signal  Gradient 

%  1  volt/degree 

Input  Saturation  Level 

%  2.5  V  RMS 

Output  Level 

2  Amps  dc  (max),  300  V 
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4.5.2  Acoustic  Parametera 

Transducer  Material 
Diameter 
Directivity 
Beamwidth  (3  dB) 
Side  Lobes  (first) 
Efficiency 
Power  Output 
Receiving  Sensitivity 


Barium  Titanate  Disk 

4  inches 

*^30  dB  actual 

6  degrees 

18  -  22  dB  down 

~  50% 

1  watt/projector 

~  -83  dB  vs  1  volt/pbar  at  1  yard 
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5.0  CONCLUSIONS 

The  research  and  experimental  verification  thereof  have  proven  that  a  suitably  con¬ 
figured  Doppler  navigation  system  can  provide  extremely  accurate,  bottom-referenced, 
speed  and  position  information.  In  the  course  of  this  study,  certain  other  areas  of  un¬ 
certainty  have  become  evident  although  they  in  no  way  should  cast  doubt  upon  the  validity 
of  the  work  already  accomplished.  These  areas  are;  1)  The  need  for  accurately  deter¬ 
mined  volume  reverberati(m  strengths  and  2)  A  more  accurate  description  of  the  back- 
scattering  process.  References  1,  6  and  7  discuss  these  problems. 
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